Galaxy cluster number count data constraints on dark energy 
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Abstract 

We use data on massive galaxy clusters (Mdustcr > 8 x W^^h~^ Mq within a comoving radius of 
J^ciuster = 1.5/i~^Mpc) in the redshift range 0.05 < z < 0.83 to place constraints on the dark energy 
equation-of-state parameters (wojWa), as well as on the nonrelativistic matter density parameter 
Q,m., on the amplitude of mass fluctuations erg, on the index n of the power-law spectrum of the 
density perturbations, and on the Hubble constant Hq. We show that present data weakly constrain 
{wo,Wa) around the values corresponding to a cosmological constant, i.e. {wo,u)a) ~ (—1,0). Clus- 
ter data alone prefer low values of the amplitude of mass fluctuations, as < 0.69 (Icr C.L.), and large 
amounts of nonrelativistic matter, Qrn > 0.38 (Icr C.L.), in slight tension with the ACDM concor- 
dance cosmological model (the results are however compatible with ACDM at 2a). We derive a erg 

1/3 

normalization relation, crgflm ~ 0.49 ±0.06 {2a C.L.), and we find that the fit to data is almost in- 
dependent of n and only very weakly dependent on Hq. Combining cluster data with baryon acoustic 
oscillation observations, cosmic microwave background data, Hubble constant measurements, Hub- 
ble parameter determination from passively-evolving red galaxies, and magnitude-redshift data of 
type la supernovae, we find erg = 0.731^^^, flm = 0.28tH2> = -l.Uto.lt^ = 0.85tofo, 



Ho = 69.1+i j km/s/Mpc (all la C.L. errors), in substantial agreement with the concordance cos- 
mological model. Very similar results are found in the case of time-evolving dark energy with a 
constant equation-of-state parameter w — const (the XCDM parametrization). Finally, we show 
that the impact of bounds on {wo,Wa) is to favor top-down phantom models of evolving dark energy. 
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PACS numbers: 95.36.-l-x, 98.80.Es 



I. INTRODUCTION 

In the last few years galaxy cluster observations have 
begun to provide useful constraints on cosmological pa- 
rameters (for a recent review see Ref. [H). Data on galaxy 
clusters are now used to test the validity of the stan- 
dard cosmological model, the so-called A cold dark mat- 
ter (ACDM) concordance model, [2|, which describes ob- 
servational data at large cosmological scales (from galac- 
tic scales to the present horizon scale) reasonably well 
0. In particular, cluster observations can help tighten 
the bounds on cosmological parameters such as the non- 
relativistic matter density parameter f2„i, the amplitude 
of mass fluctuations ag, the power-law index n of the 
density perturbation power spectrum, and the Hubble 
constant Hq Q. 

When combined with other cosmological probes — 
such as cosmic microwave background (CMB) radiation 
anisotropy, baryon acoustic oscillations (BAO) in the 
matter power spectrum, Hubble parameter, and type la 
supernovae (SNela) data — galaxy cluster observations 
provide a unique insight towards helping understand the 
evolution of the Universe, from the inflation era to today. 

Despite the observational success of the ACDM model, 



a number of basic questions remain unanswered. Dark 
energy is a major mystery (for reviews on dark energy 
and modified gravity see, e.g., Refs. [1]). A possibility is 
that dark energy is simply a manifestation of a nonzero 
vacuum energy, a cosmological constant Aj but dynam- 
ical scalar field models of dark energy, 0|, are also 
compatible with present data.^ 



It has been known for some time that a spatially-flat ACDM 
model is consistent with most observational constraints, see, e.g., 
Refs. ^6). In the ACDM model the energy budget is dominated 
by far by a cosmological constant, a spatially homogenous fluid 
with equation of state parameter w\ = pa/ PA = — 1 (where pA 
and Pa are the fluid pressure and energy density), with nonrela- 
tivistic CDM being the second largest contributor. Note that the 
"standard" CDM structure formation model — which the stan- 
dard ACDM cosmological model assumes — might have some 
observational inconsistencies (see, e.g., Q). 

In dynamical dark energy models the dark energy density de- 
creases in time and so remains comparable to the nonrelativis- 
tic matter (CDM and baryons) energy density for a longer time 
(than does a time-independent A). This partially alleviates the 
"coincidence" puzzle of the ACDM model. In addition, some 
dynamical dark energy scalar field models have a nonlinear at- 
tractor solution that generates the current, tiny, dark energy 
density energy scale of order an meV from a significantly higher 
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Measurements of the local abundance and growth of 
galaxy clusters from X-ray l3"13| and optical [31 sur- 
veys have been recently used to probe the standard cos- 
mological model. In particular, the emerging picture is 
that a cosmological constant still remains a good candi- 
date for dark energy. This conclusion does not exclude 
the possibility that future cluster surveys will allow us to 
discriminate between the ACDM model and dynamical 
dark energy models [ll] . 

The aim of this paper is to use present data on massive 
galaxy clusters (Mdustcr > 8 x IO^'^H'^Mq within a co- 
moving radius of -Rciuster — 1.5/i~^Mpc) at low and high 
redshifts (0.05 ^ z < 0.83) to constrain some of the free 
parameters of the standard cosmological model, and to 
investigate the possibility that the dark energy density 
evolves in time, instead of staying constant. In Refs. [16[ 
these data were used to determine the linear amplitude of 
mass fluctuations and the nonrelativistic matter density 
in a Universe with a cosmological constant. We extend 
the analysis of Refs. [l^ to the case of evolving dark en- 
ergy, and we properly take into account the dependence 
on redshift and cosmology of quantities related to cluster 
physics: the critical density contrast, the growth factor, 
the mass conversion factor, the virial overdensity, and the 
virial radius. 

The plan of our paper is as follows. In the next section 
we briefly discuss some models of evolving dark energy 
and introduce the parametrization of the time-dependent 
equation of state we use in our analysis. In Section III, 
after introducing the basic theory of galaxy cluster num- 
ber counts, we use data from X-ray observations of mas- 
sive clusters to place constraints on cosmological param- 
eters. In Section IV we present data from other types of 
cosmological probes: baryon acoustic oscillations, cosmic 
microwave background radiation anisotropics, passively- 
evolving red galaxies, and type la supernovae. In Section 
V we outline a joint analysis of all data, while in Sec- 
tion VI we describe our conclusions. Finally, in the Ap- 
pendices we discuss in detail the critical density contrast 
and the growth factor (Appendix A), the mass conver- 
sion factor (Appendix B), and the virial overdensity and 
the virial radius (Appendix C). 



II. DARK ENERGY MODELS 

It is of great interest to determine if the dark energy 
is well- approximated by a cosmological constant or if it 
decreases slowly in time (and so varies weakly in space). 



Ideally one would very much prefer a mo del- independent 
resolution of this issue. However, at this point in time, 
observational data are not up to this task. One must 
instead use the available observational data to constrain 
model parameters and so determine if the cosmological 
constant point in model parameter space is or is not fa- 
vored over points where the dark energy density slowly 
decreases in time. While it is useful to perform such 
an analysis using a consistent and physically motivated 
model, such as the inverse power-law potential energy 
density scalar field model d, 0, this is computation- 
ally quite demanding, so here we make use of a simple 
parametrization of time-evolving dark energy in a pre- 
liminary attempt to investigate this matter. 

In order to discriminate between a cosmological con- 
stant and dynamical dark energy we use the dark energy 
equation-of-state parameter parametrization [l7| 



w{a) = wq + Wa {I — a) 



(1) 



where a is the scale factor related to the redshift z by 
a = 1/(1 -I- z). The cosmological constant corresponds 
to Wo = —1 and Wa = 0, the case of constant equation 
of state corresponds to wq — w — const and Wa — 
(known as the XCDM parametrization of time-evolving 
dark energy) , while the general case of time-evolving dark 
energy corresponds to Wa ^ 0. 

The advantage of the above parametrization is twofold: 
(i) a number of dark energy models can be adequately 
described by an equation of state of the form ([T]) at recent 
enough times (i.e., for a near unity); and, (ii) at a given 
redshift (such that a ^ 1) different classes of dark energy 
models correspond to different regions in the {wo,Wa) 
plane, [3]. Indeed, roughly speaking, there exist four 
classes of dark energy models: "thawing" models, "cool- 
ing" models, "barotropic" fluids (all assumed to obey the 
null energy condition w > —1), and "phantom" models 
(for which w < —1).^ Introducing the quantity 



w 



dw 
din Q 



(2) 



the classification is as follows. 

Thawing models - These satisfy the inequalities [l9| 

l + w <w' <3il + w), (3) 

and can arise in models of dark energy implemented by 
(cosmic) scalar fields, such as axions or dilatons, which 
roll down towards the minimum of their potential. Typ- 
ical potentials are of the form (p™ (to > 0), with (j> being 
the scalar field. The bounds ^ are valid for (1 + w) <C 1 



energy density scale (possibly of a more fundamental model) as a 
consequence of the very slow decrease in time of the dark energy 
density during the very long age of the Universe. These results 
are often viewed as providing significant theoretical motivation 
to consider dynamical dark energy models, 



^ This classification is not exhaustive since, as explained in llSlI , 
some models, such as those with oscillating equations of state, 
do not fall into it. 
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so, following Ref. [l9[, we assume it; < —0.8 as a practical 
limit of applicability. 

Cooling models- As for the case of thawing models, 
cooling models may arise in scalar field models of dark 
energy. Typical scalar potentials are of the form <j>~™ 
(m > 0). They lie in the region [11, 



-3{l-w){l + w)<w'<e{z)w{l + w), (4) 

where e is a function of the redshift and an investigation 
of a variety of scalar-field cooling models indicates that 
e(l) ~ 0.2 The upper bound in Eq. (jU is vahd for 
{1 + w) <^ 1 [l^ so, for this bound, we assume w < 
—0.8. Cooling models may arise in models of dynamical 
supersymmetry breaking and supergravity. So-called k- 
essence models [2l| with a nonlinear kinetic term belong 
to this class [2^ . 

Barotropic fluids - These are fiuids whose pressure p 
depends only on the energy density p. Assuming that 
Cg = dp /dp > 0, barotropic fiuids satisfy the inequal- 
ity 



w' < 3u;(l -I- w) , 



(5) 



and include (original [23| and generalized [24| ) Chaplygin 
gas models. 

Phantom models.- These are models which do not obey 
the null energy condition (see, however, Ref. 18|), and 
fall into the region 



w < -1 



(6) 



To each of the above models, we can associate a spe- 
cific region in the (wQ.Wg) plane at a given reference time. 
Following Ref. |l8| , we take as reference time that corre- 
sponding to z = 1 which, roughly speaking, is when dark 
energy is expected to start to dominate over nonrelativis- 
tic matter. Phantom models at z = 1 can be split in two 
classes: "Pure phantom" models which did not cross the 
phantom divide line w ~ —1 recently. 



wq < —1, w{z = 1) < — 1 



(7) 



and models that crossed w = — 1 from a lower value to a 
higher value. 



Wo > —1, w{z = 1) < — 1 



(8) 



The {wQ,Wa) plane containing all the above regions is 
presented in Fig. 4, together with the regions allowed by 
data and discussed in Section V.'' 



III. GALAXY CLUSTER NUMBER COUNTS 

In this section we introduce the basic physical quan- 
tities and observables related to galaxy cluster number 
counts, we analyze the available experimental data, and 
we discuss the results of our analysis. 



Ilia. Theory 

Cluster number and comoving volume - The comoving 
number of clusters in the redshift interval [zi,Z2], whose 
mass M is greater than a fiducial mass Mg, is 



'dz^^^N(M > Mo,z) 
dz 



where 



V(z) = An fdz' "^'jf^ , 

^' Jo il + z'^Hiz') 

is the comoving volume at redshift z, and 

dz' 



dUz) = (1 + z) 



H{z') 



(10) 



(11) 



(12) 



is the luminosity distance with H(z) the Hubble param- 
eter. The "mass function" N(M > Mo,z) appearing in 
Eq. (jlOp represents the comoving cluster number density 
at redshift z of clusters with masses greater than Mq. 

For a cosmological model with evolving dark energy 
equation-of-state parameter of the form ([T]), the Hubble 
parameter normalized to its present value Hq is 



E{z) = 



H{z) 



Per 



Per 



1/2 



The quantities 



p„,{z) = f7„p(?(l + z) 



(13) 



(14) 



which we dub "bottom-up phantom" models. 

Finally, we also consider models that crossed w = —1 
from a higher value to a lower value: 



wq < —1, w{z = 1) > — 1 



(9) 



These models, which we dub "top-down phantom" mod- 
els, are phantom today (z = 0) and non-phantom at 
z=l. 



* There are a few disadvantages of the parametrization of Eq. 
|[TJ: (i) This two parameter (wo,Wa) parametrization has one 
more parameter than the simplest consistent and physically mo- 
tivated scalar field dark energy model, [^l9j, thus making it more 
difficult to constrain model parameters with observational data; 
and, (ii) even so, the parametrization is not physically complete 
as additional information must be provided if one is interested in 
the evolution of energy density and other spatial inhomgeneities, 

M. 
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and 

PDe(^) = f^DEpi?(l + z)3(l+-o+-.) g-3«,„./(l+.) (15) 

are the energy densities of nonrelativistic matter and 
dark energy, respectively. Here, Qm = Pm / Pcx and 
f^DE = Pbe/Pct:'' matter and dark energy den- 

sity parameters, and pm\ Pqej and pc?'' = 3i?g/(87rG) 
are the present matter, dark energy, and critical energy 
densities, respectively, while G is the Newton constant. 

In this paper, for computational simplicity, we restrict 
ourselves to the case of a flat universe,^ so 



(16) 



Mass function - To compute the mass function, we use 
the Press-Schechter (PS) ap proach ^28.] , as modified by 
Sheth and Tormcn (ST) [29|. In this approach, the mass 
function is written as 



N{M>Mo,z)= dMn{M,z), (17) 

Jmo 



and 



M 



47r 



'r^{z)Pm{z) 



An 



(18) 



is the mass within a sphere of physical radius r{z), whose 
corresponding comoving radius is i? = (1 + z) r{z). 

In Eq. (jl7p . n(M, z) dM is the comoving number den- 
sity at redshift z of clusters with masses in the interval 
[M,M -I- dM], and is written as 



2/Om , dv 

-W^^'^dM 



(19) 



Here the multiplicity function i^/(i^) is (in the PS and ST 
models) an universal function of the peak height 



Sc 

V = — , 



and is normalized as 



r°° 1 

dvvf{v) = - 



(20) 



(21) 



The functional form of vfiv) is discussed below. The 
critical density contrast 5c{z) is the density contrast for 
a linear overdensity able to collapse at the redshift z, and 
its dependence on cosmological parameters is discussed 
in Appendix A. 

The root mean square (rms) amplitude a of density 
fluctuations in a sphere of comoving radius i?, whose 



This is consistent with the simplest interpretation of the CMB 
anisotropy data, see, e.g., [26l l27ll . 



corresponding physical radius r contains the mass M, 
is related to the matter power spectrum of density per- 
turbations at redshift z, P{k,z), through 

1 r°° 

a^{R,z) = — dkPP{k,z)W^{kR) . (22) 
271-^ Jo 



Here 



W{x) = 



3 (sin a; — xcosx) 



(23) 



is the Fourier transform of the top-hat window function 
and 



P{k,z) = Po{k)T^{k)D'^{2 



(24) 



with D{z) being the growth factor (discussed in Ap- 
pendix A) and T{k) the transfer function. 

We assume that the post-inflationary density pertur- 
bation power spectrum Po{k) is a simple power law. 



Po{k) = Ak'' 



(25) 



with the scale-invariant spectrum corresponding to n = 
1. The normalization constant A is a free parameter of 
the model and can be expressed as a function of the other 
cosmological parameters (see below, footnote 6), while 
the total transfer function T(k) is taken from Ref. f3Q |. 
The transfer function depends on Hq and on baryon and 
cold dark matter density parameters flf, and Q^- The 
total amount of matter is given by fim — Qb + and in 
this paper we take Qi,h^ = 0.02, with h defined by 



Ho = lOOh km/s/Mpc 



(26) 



In the Press-Schechter parametrization, the form of 
the multiplicity function is a result of the assumption 
that initial density fluctuations are Gaussian: i'f{i') = 



/ \j2/K. In this paper, however, we use the form 
v!(v)^K{\ + {av')-^\e-'^^'l\ (27) 

introduced by Sheth and Tormen, inspired by a model of 
elliptical collapse. The constant 



K = (28) 

V2^+2i/2-Pr(l/2-n) 



2i/2-pr(i/2-p) 

results from the normalization condition (1211) . r(a;) is the 
Gamma function, while a and p are phenomenological 
constants to be determined by fitting to A'^-body simu- 
lation results. We use the values found by Sheth and 
Tormen, namely a = 0.707 and p = 0.3. (The Press- 
Schechter case is recovered for a = 1 and p = 0.) 

Finally, putting all this together, we can rewrite the 
mass function as 



A(M > Mo,z) = 

Pmiz) Sc{z) 



(29) 



2K 



Mn 



c7sD{z 
aS^iz) 



•dx 1 



Y?{x) 



dEix) 



dx 



cr|L>2(z)E2(2.) 









> exp 



aSliz) 



2alD^{z)j:^{x) 
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where we have introduced the function 

. _ /r dyy"+^ T\y/Rs) W^xyRp/Rs) 
^""^ S^dyy-+^T^{y/R,)W^{y) ' ^'"^ 

and the quantities^ 



(78 = cr(^8,0) , 
i?8 = 8/i"^Mpc 



(31) 



3A/o 



1/3 



(0) 



We note that the function evaluated at x = R/Rq 
is the present value of the rms amplitude a at the scale 
R normalized to its present value at the scale Rs: 



S(i?/i?o) = 



a{R,0) 



^8 



(32) 



This result will be used in Appendix A. 



Illb. Data 

Data on cluster abundance at different redshifts can be 
expressed as the comoving number of clusters in the red- 
shift interval [zi , 22] , with mass M' within a reference co- 
moving radius R'q greater than a fiducial mass Mq within 
the same radius. Here and in the following, a prime is 
used to distinguish quantities related to observed masses 
and radii from theoretical ones, discussed in the previous 
subsection. We follow Refs. 161 and take 



R'q = 1.5/i"^Mpc , 
Afo = 8 X 10"/i~^Mp 



; 



(33) 
(34) 



where Mq ~ 1.989 x 10 g is the solar mass. 

Only an effective fraction a{z) of the total comoving 
volume at redshift z is observed, so the expected comov- 
ing number of clusters in the redshift interval [21,^2], 
with mass M' greater than Mg, is 



Jz, dz 



where N'{M' > Mq, z) represents the comoving cluster 
number density at redshift z of clusters with masses M' 
greater than M(). The mass function iV'(M' > M[t,z) 
can be written as 

/•OO 

N'{M'>Mq,z)^ dM'n'{M\z) . (36) 



As anticipated, the normalization constant A can be related 
to the other cosmological parameters, f!;,, f2c, and n: A = 
2n^4R-+^/ dyy"+^ T\y/Rs) W^y). 



TABLE I: The four redshift intervals (« = 1. 2, 3, 4), 

centered at z'*' , of the massive clusters data and the references 
from which data have been taken, ai is the effective fraction 
of the observed comoving volume of the i*'' bin and the values 
listed here are for the case with = 0.3 and w = —1. 



bin i 




^1 




Ref. o-i 


1 


0.00 


0.10 


0.050 


[31] 0.309 


2 


0.30 


0.50 


0.375 


[32] 0.012 


3 


0.50 


0.65 


0.550 


[16] 0.006 


4 


0.65 


0.90 


0.825 


[33] 0.001 



Here n'(M' ,z)dM' is the comoving number density at 
redshift z of clusters with masses M' in the interval 
[M', M' H- dM% and is defined by 



n'(M', z) dM' = ?i(M, z) dM , 

where n(M, z) is given by Eq. (fT9|) . 

Inserting Eq. ([57)1 in Eq. (1551) we obtain 



(37) 



N' [W > M'q, z) ^ I dMn{M,z), (38) 

where the function g relates the observed mass M' to the 
virial mass M in the PS or ST parametrization. Conse- 
quently, g^M^) is the fiducial virial mass 



Mo = g(M^) 



(39) 



which corresponds to the fiducial mass Mq adopted in the 
observations. In Appendix B we describe the procedure 
that gives the mass Mq as a function of Mq. In general, 
the function g depends on the redshift and cosmological 
parameters. 



In Table I, we show the four redshift bins 

At) 



{i = 1,2,3,4), centered at Zc' , of the massive clusters 
data. Also listed are the values of the effective fraction 
of the observed comoving volume of each bin, a^. The 
values are for a cosmology with {VL„n w) = (0.3, —1), and 
were computed using the results of Refs. |16{ . 

The ai parameters depend, in principle, on the cosmol- 
ogy and their values can be obtained using the E(l / Knax) 



method applied to observational data [31|, |32|, |3J| . How- 
ever, the dependence of ai on the cosmology is weak com- 
pared to that of the comoving volume. Using the results 
of Refs. 16[ we get, for example, that passing from the 
cosmology with (fi^, w) = (0.3, —1) to that with = 1, 
the percentage variation of the comoving volume relative 
to the third bin, T^(4^^) - Vizi^), is 117%, while that of 
as is about 1%. Similar results hold for the other bins. 
References 

0, 

[3ll - l33j give X-ray temperature mea- 
surements for massive clusters. For completeness we 
show these data in Table II. 

In order to convert temperature to mass, we use the 
mass-temperature conversion formula of Ref. (ssj (see 
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TABLE II: Name and X-ray temperature Tx of clusters in the 
four bins used in our analysis. Data in the first and second 
bins are from Ref. [3l[ and Ref. [S^l, respectively, while data 
for the third and fourth bins are from Ref. [l^ and Ref. [s^l, 
respectively. All errors are at the 68% confidence level. 



bin i 


name 


Tx (keV) 




A0754 


9.oot;5;li 




A2142 






COMA 






A2029 






A3266 


_ „r, + 0.21 




A0401 


— ^ r\-\-Q 17 

7.i9Io;;5 




A0478 






A2256 


c OO+0.14 




A3571 






A0085 


„ - _i_n in 




A0399 


6.46^^:^^ 




ZwC11215 


/-> 0/^ + 1.79 

6.36Iii^ 


1 


A3667 


/-> oo4-0 16 




A1651 


6.22l°:^^ 




A 1795 


^ - rj-4-0 16 




A2255 


c nn+0.24 




A3391 


c QQ+O 27 
t).89_0 20 




A2244 


0- ' '-0.27 




A0119 


5.69^^.1? 




A1650 


r co + O.lS 
0.0»_o.i6 




A3395s 


r rc;+0.54 
0.00_o 40 




A3158 


f- ,.,+0.16 

5-41Io.i5 




A2065 


r n>7+0.21 
0-'J'-0.18 




A Q C CO 


r 07+0.10 
^-"J'-0.09 




A3112 


4-72Io.i5 




A1644 


4 70+n,^ 
^- ' ^—0.43 




MS 1008.1 


8.2II:? 




MS 1358.4 




2 


MS 1621.5 


"-"-0.8 




MS 0353.6 


6-5lJ:^ 




MS 1426.4 


6.411:° 




MS 1147.3 


e-otJ:? 


3 


MS 0451-03 


10.41°:^ 




MS 0016+16 


0+0.6 

8-0.6 


4 


MS 1054-03 


12.311:^3 



Ref. [36[ for a different approach to the problem of cluster 
mass-temperature conversion): 



Af (< r) = W^Mq ka 



Tx 



keV Mpc 



(40) 



where M(< r) is the mass within a physical radius r, Tx 
is the cluster X-ray temperature, and ka is a parameter 
which depends only on A^. Here is the virial over- 
density relative to the critical density. It is related to 



the virial overdensity relative to the background matter 
density, A.„, through 



(41) 



The quantity A^, depends on redshift and cosmology and 
is thoroughly discussed in Appendix C.^ 

As found in Ref. [35|, the parameter ka depends on 
A^ and, in particular, when 

a; = 100, 250, 500, 1000, 2500 , (42) 

ka assumes, respectively, the values 

KA = 0.76, 0.91, 1.01, 1.09, 1.14 . (43) 

From Eq. (|40p . it follows that a mass M' within a co- 
moving radius Rq = 1.5/i~^Mpc is related to the X-ray 
temperature by 



M' = 1.5 X lO"/i"^M0 KA 7-^ 



Tx 1 



keV 1 + z 



(44) 



This means that clusters in the i"^ bin, with masses 
Af > A/q, will have a temperature exceeding the thresh- 
old value 



Tx,o 16(1 + z'c 



keV 



3ka 



(45) 



In order to apply the above equation to data, we ex- 
trapolate the parameter ka according to the following 
prescription: 



0.76 if A'„ e [25, 175] , 

0.91 if A'„ G ]175,375] , 

1.01 if a; g ]375,750] , 

1.09 if a; e ]750, 1750] , 

1.14 if a; e ]1750,3250] 



(46) 



By using Eqs. (|45|) and (|46| and data from Table II we 
find the values listed in Table III for the observed number 
of clusters with masses M' > AIq in the i*^ bin, A/'^j^g j. 

The uncertainty in the comoving numbers of clusters, 
AA^Qj^g j , derive from the uncertainty in the X-ray tem- 
perature of clusters. The threshold X-ray temperature 
in each bin and for each A^ interval is also indicated in 
Table III. 



Defining as the physical radius containing an overdensity 
of relative to the critical density, and M{r < r^) as the 
mass contained in rA, the mass-temperature relation II40I I as- 
sumes the standard T^''^ power-law form [35l|, namely M{r < 

[3/4.A;p(?l^/^l + .)-3/2Tr. 
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TABLE III: The threshold X-ray temperature, Tx.o, and the observed number of clusters with masses M' > Mq, M^hs ii in 
each bin and for each interval. The uncertainty on the comoving numbers of clusters, AA/'obs.ii is also indicated. 







bin 1 


bin 2 


bin 3 


bin 4 






[rx,o(keV),Kb,,i] 


[Tx,o(keV),Xbs,2] 


[Tx,o(keV),Xb,,3] 


[rx.o(keV),AA^, 


a; 


£ [25, 175] 


[7.37, 5+J] 


[9.6, Ol°] 


[10.9, OlJ] 


[12.8, OlJ] 


a; 


€ ] 175, 375] 


[6.15,151^] 


[8.1 , 11?] 


[9.1, llj] 


[10.7, ll°] 


a; 


G ]375, 750] 


[5.54,211^5] 


[7.3,111] 


[8.2,21?] 


[9.6,11^] 


a;, 


e ]750, 1750] 


[5.14,2411] 


[6.7, 2lt] 


[7.6,21?] 


[8.9, ll°] 


A^ 


G ]1750, 3250] 


[4.91, 24lg] 


[6.4, 5lJ] 


[7.3, 2l°] 


[8.5, ll°] 



Finally, we calculate the observed number of clusters 
in the four redshift bins as 



dz^^N'iM' > M^,z) , (47) 



where we used Eq. ([55)1 and replaced a{z) with the aver- 
age effective fraction ai in the i"^ bin. Since, as argued 
above, the ai only depend weakly on the cosmology, we 
use in Eq. (H71) . for definiteness, the values for the case 
rim — 0.3 and w — —1 listed in Table I. 

In Fig. 1 wc plot the various quantities needed in the 
computation of the mass function as a function of redshift 
and for different choices of {wq, Wa)- We fixed the values 
of the other cosmological parameters to the best-fit values 
obtained by using the 7-year WMAP observations '27^, 
namely {h,n,n„,,as) = (0.71,0.96,0.3,0.80). Although 
the comoving volume is very sensitive to the choice of 
the cosmological model, the variations of the functions 
6c, D{z), Mq/Mq, and A„ (discussed in the Appendices) 
with redshift and cosmology concur to give rise to larger 
changes in N'{M' > Mq, z), especially at large redshifts. 
This, in principle, can be used to put constraints on var- 
ious models of dark energy. 



IIIc. Results 

Due to the small number of clusters in each bin, the 
comparison between observed and predicted number of 
clusters is made using Poisson error statistics. Accord- 
ingly, we define a likelihood function by 



(48) 



where we have introduced = A/"/ and Ki = A/'^i^g ^ for 
notational clarity. The statistics is then introduced as 

XCL,nosys(^'"-'^m'O-8,W'0,Wa) = -2lnC (49) 
4 

^2J2[K~ K^.s,^ (1 + In A/^' - In A/^^s,.)] ■ 



We also take into account the uncertainty in the comov- 
ing numbers of clusters, AA'^i-.g j, by shifting the observed 
number of clusters in each bin as 



obs 



obs,^ 



eAAA.^bs.^^A/'o'bs 



(50) 



where the "pull" ^ is a univariate gaussian random vari- 
able [s^]. Correspondingly, we modify the as 

4 

= 2 ^ [A/;' - AA^bs,. (1 + In A/;' - hi AA^b,,,)] + e ■ 

i=l 

Since the depends on seven parameters, a grid-based 
analysis is not feasible and we therefore employ a Markov 
Chain Monte Carlo approach. We use a modified version 
of CosmoMC ^] to produce and analyze the likelihood 
chains. 

Figure 2 shows the results of the analysis in the 
{wo,Wa) and {flrn,<^s) planes. Here, we marginalize over 
the other parameters, using a fiat prior, and determine 
the Icr, 2a, and 3cr regions by finding where increases 
by 1, 4, and 9, respectively, from the most likely set of 
parameters. 

The filled contours are obtained taking into account 
only statistical uncertainties (i.e., taking for Af^y^^ ^ just 
the best fit values listed in Table III), while empty con- 
tours show the effect of including systematic errors on 
the comoving numbers of clusters. 

The poor quality of current data on massive clusters 
do not allow us to appreciably constrain the equation-of- 
state parameters {wQ,Wa), to either favor or rule out a 
cosmological constant as dark energy. For this reason, in 
the next section we combine cluster data with data from 
other cosmological probes to better constrain {wQ,Wa)- 

We find the marginalized bounds (including systemat- 
ics), 



n,n > 0.38 (Icr C.L. 
as < 0.69 (Icr C.L. 



(52) 



These results can be compared to other recent cosmolog- 
ical results from galaxy clusters obtained using different 
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FIG. 1: The comoving volume V{z), the critical density contrast 5c{z), the growth factor D(z), the mass ratio Mq/Mq, the 
mass function N'{M' > Mq,z) [see Eq. ((38])], and the virial overdensity relative to the background matter density, A„, as a 
function of the redshift z, for different dark energy equation of state parameters (wo,iOa). We fixed the values of the other 
cosmological parameters to (/i, n, fim, fg) = (0.71,0.96,0.3,0.80). 



strategies and cluster surveys. Studies of X-ray selected 
clusters, with masses exceeding a fixed mass threshold 
and distributed over fixed redshift ranges, yield lower val- 
ues of rim and larger values of cts compared to Eqs. ((52)) 
in better agreement with the concordance cosmological 



model emerging from data analysis of other cosmologi- 
cal probes, such as CMB anisotropy and SNela [ij. For 
example, Mantz et al. [HI find Vtm = 0.23 ± 0.04 and 
CTs = 0.82 ± 0.05 in a model with a constant dark energy 
equation of state w. However, both the analysis of cluster 
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FIG. 2: la, 2a, and 3a confidence level contours in the {wQ,Wa) and {Q,m,as) planes, determined from galaxy cluster num- 
ber count data. Results obtained by including systematic uncertainties are shown as empty contours, while those ignoring 
systematics (i.e., keeping just best fit values for data) are represented as filled contours. 



population over low redshift ranges [39[ and the study of 
the evolution of the Tx-based mass function [l^] , give re- 



sults compatible with ours, namely flm = 0.34 



as = 0.71 



-1-0.13 
-0.16' 



-1-0.09 
-0.08 



and 



and flm = 0.34±0.08, respectively (these 



results are obtained assuming a model with a cosmolog- 
ical constant). 

Our fit to data gives the erg normalization 



1/3 



0.49 ±0.06 (2(tC.L. 



(53) 



This should be compared to the results of Refs. [16[ ob- 
tained from an analysis of the same data we use in this 
paper:8 asfl^J = 0.33 ± 0.03 and dsf^?™" = 0.78 ± 0.08, 
resulting from fit to data relative only to the first and last 
two bins in Table I, respectively. Our 2cr-band ([55)1 in the 
(Tg — flm plane lies between the two bands defined by the 
above two relations for flm ^ 0.4, with the three bands 
overlapping in the region near flm — 0.2. So, we con- 
clude that our ag normalization agrees with the results 



In Refs. 16], however, it is assumed that h = 0.72, n = 1, and 
w = —I. Moreover, the parameter in Eq. II40I I is assumed 
to be cosmology-independent and the value used is simply the 
arithmetic mean of the values in Eq. II43I I. namely ka = 0.98. 
Also, the mass conversion is done by using the observed cluster 
profile in the comoving radius range R G [0.5, 2] h^'^Mpc, while 
we assume a Navarro-Frenk- White profile for the virialized halo 
mass density (see Appendix B). Finally, the expressions for the 
critical density contrast, growth factor, virial overdensity, and 
virial radius refer, in Refs. [16| |. to a matter dominated universe 
with Qm = 1. 



of Refs. [16 1 . 

We also find that the results of the fit are almost in- 
dependent of n in the adopted range [0.90,1.05], and 
that they are only very weakly dependent on h in the 
range [0.6,0.8]. Finally, the pull ^ is found to be ^ 
-0.2ll°:i^ (laC.L.). 



IV. OTHER COSMOLOGICAL DATA 

In this section, we present data from other type of 
cosmological observations. In the next section, we de- 
rive joint constraints using these data along with those 
of massive clusters. 

Baryon Acoustic Oscillations - The measurement of 
baryon acoustic oscillations (BAOs) in the large-scale 
matter correlation function, fixes the values of a char- 
acteristic "BAO distance parameter", which we denote 
asC.9 

With Dv an effective distance defined by 



Dv{z) 



1 



( z V^^ 


■ r dz' ' 


2/3 




[Jo E{z')\ 





(54) 



® See Refs. [4Cll | for recent discussions of BAO data constraints on 
cosmological parameters. 
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the C parameter is the ratio 



C 



Dv{0.275) 



(55) 



between the comoving sound horizon at the baryon drag 
epoch Zd, 



Tsiz) 



1 



dz' 



417^(1 + z) 



ViHoJz E{z')\ An^{i + z) + 3nb 



(56) 



and the effective distance Dy at z = 0.275 [4l|. Here, 
flj is the photon density parameter that we take equal 
to n-yh'^ = 2.56 X 10"^ [H]. 1° For the redshift at the 
baryon drag epoch, Zd, we use the fitting formula of Ref. 



Zd = 



where 



1291(f7„/i2)C 



[l + 6i(r!,/^2)h. 



1 + 0.659(r!™/l2)0.828 



(57) 



h - 0.313(an/i')"°-^'^ [1 + 0.607(r!„,/i2)°-^^^] , (58) 
62 = 0.238(17„,/i2)0-223 (59) 

BAO data give the value Cobs = 0.1390 ± 0.0037 
Accordingly, we define the statistic 

2 (C - 0.1390)^ 

XbAoI^i ^m,Wo,Wa) = ^ ^^^^^ ' (^O) 

Cosmic Microwave Background - The analysis of the 
CMB radiation puts a constraint on the reduced distance 
to the surface of last scattering, the so-called "CMB shift 
parameter" , 



' dz 

W) 



(61) 



where zis ~ 1090 is the redshift at the time of last 
scattering. The shift parameter weakly depends on the 
adopted cosmology and here we use the constraint found 
by Corasaniti and Melchiorri, 7^obs = 1-710 ± 0.026 Q, 
which refers to a cosmological model with evolving dark 
energy with equation-of-state parameter of the form 
given in Eq. ([T]). We then consider the following statis- 
tic 



XCMBi^m, Wo, Wa) 



(7e- 1.710)^ 
0.0262 



(62) 



Since the upper limit of integration of the integral in Eq. II56I I 
is infinity we must include, in the expression of the normalized 
Hubble parameter E{z), the contribution due to radiation. Ac- 



cordingly, we add the quantity pr / p, 



(0) 



: nr(l + z) in the argu- 



ment of the square root appearing in Eq. I I13II , where pr and £7,. 
are the radiation energy density and radiation density parameter, 
respectively We take Qrh^ = 4.31 X IQ-^ [il . 



TABLE IV: The observed Hubble parameter Hobs{zi) with 
error an (in brackets) from passively evolving galaxies (data 
from Ref. and line-of-sight BAO peak position (data are 
from Ref. [49| and marked with an asterisk). 



Zi 


Hobsizi) [km/s/Mpc] 


0.1 


69(12) 


0.17 


83(8) 


0.24 


79.69(2.65)* 


0.27 


77(14) 


0.4 


95(17) 


0.43 


86.45(3.68)* 


0.48 


97(60) 


0.88 


90(40) 


0.9 


117(23) 


1.3 


168(17) 


1.43 


177(18) 


1.53 


140(14) 


1.75 


202(40) 



Hubble Constant - A meta-analysis of many measure- 
ments yields Hq = (68±2.8) km/s/Mpc at la C.L. [H." 
Accordingly, we introduce the penalty 

Xlih) - "'-TJ' . (63) 



(fe- 0.68) ^ 
0.0282 

Hubble Parameter.- The analysis of spectra of passive- 
ly-evolving red galaxies enables the determination of the 
Hubble parameter at different redshifts, [i^]. We use 
data quoted in Ref. [481 and reported in Table IV for 
the sake of completeness. To these data we also add the 
estimate of the Hubble parameter at redshifts z — 0.24 
and z = 0.43, obtained in Ref. (49j by using the BAO 
peak position as a standard ruler in the radial direction. ^2 

We then introduce a statistic as 
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[H{zi) - Hohsizi)] 



(64) 



'H 



Type la supernovae - Type la supernovae are standard- 
izable candles and so can be used to discriminate between 
different cosmological models. Indeed, the theoretically- 
predicted distance modulus fi, defined by 



^(z) = 5 logj 



/ di 



' VlMpc 



25 



(65) 



depends on the redshift and on the set of cosmological 
parameters {h, flm, wo, Wa) and can be compared to the 



This is reasonably consistent with both 'low' [45l | and 'high' [4^ 
recent estimates of the Hubble constant. 

See Refs. [50l | for Hubble parameter measurement constraints on 
cosmological parameters. 
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FIG. 3: la, 2a, and 3a confidence level contours in the (wo, Wa) plane (left panel) and in the {Qm, as) plane (right panel) deter- 
mined from galaxy cluster number count (including systematic errors), BAO, CMB, Hubble parameter, and SNe observations. 
Filled contours are the result of combining BAO, CMB, Hubble parameter, and SN data, empty (open) thin contours give the 
confidence level contours for cluster (CL) data only, while empty (closed) thick contours are from the combination of all data. 



one "derived" from the observation of SN lightcurves [5l| , 
namely fiB- This, in turn, is deduced from the analysis 
of SN lightcurves which, if performed using the "SALT2" 
fitter 



:53, gives 
Mb 



51| 



M ■ 



(66) 



Here, m^^^ and c are the peak bolometric apparent mag- 
nitude and the color correction, respectively, while xi is a 
SALT2 fitter parameter [s^l • The absolute magnitude of 
SNe, M, and a and /? are, instead, nuisance parameters 
to be determined, simultaneously with the cosmological 
parameters from fits to data. In t his p aper, we use data 
from the Union2 SN compilation 5l| which consists of 
557 SNe. 

However, since the covariance matrices resulting from 
the lightcurve fit are not publicly available, we do not 
have any information on the correlation between the er- 
rors on rrig^^, a^i, and c. Consequently, we follow the 
analysis of Ref. 3a, 43 1 as explained in Ref. fs^ and in- 
troduce the statistic 



■Wo,Wa) 
(Mi 



A.)(^.^'-M,,)(/ifP-M,) . (67) 



The double sum runs over the 557 SNe, is the ex- 
perimental value of the distance modulus of the i*^ su- 
pernova, and 



/ii = 5 logio dL + 25 



(68) 



is the "reduced" theoretical distance modulus. The "re- 
duced" luminosity distance is 



(69) 



and afj is the covariance matrix (containing both sta- 
tistical and systematic errors), while the matrix Mij is 
given by 



E-1 - 
kl ^ik '^Ij 



(70) 



It is worth noting that is independent of the Hubble 
parameter Hq, so the in Eq. (I67|) depends only on the 
cosmological parameters {flm,WQ,Wa). 

There are many other data sets that can be used to 
constrain cosmological parameters, for example, strong 
gravitational lensing observations 55|; however, for our 
illustrative purposes here the data described above suf- 
fice. 



V. COMBINED DATA ANALYSIS 

In this section, we present the results of a joint analy- 
sis of massive cluster evolution, BAO peak length, CMB 
anisotropy, Hubble parameter, and SNe apparent magni- 
tude data. The statistic is 

(71) 



X^{h,n,njn,crs,Wo,Wa,0 



— XCL 



2 

Xbao 



2 

" XCMB 



-xl 



' XHubblc ■ 
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Wo 

FIG. 4: la, 2a, and 3a confidence level contours in the planes (woyWa) from the joint analysis of galaxy cluster number count 
(with systematics) , BAO, CMB, Hubble parameter, and SNe (the same as in Fig. 3). The shaded areas represent different 
types of evolving dark energy models according to the classification of Ref . [3 • 



and depends on the six cosmological parameters 
(h, n, fljn, as,wo,Wa), and on the pull ^. 

Figures 3 shows the result in the {wq, Wa) and erg) 
planes. Thick empty contours refer to the joint analysis, 
thin empty contours to clusters only, and filled contours 
to all data except clusters. 

The allowed la confidence limits for (wo, Wa, f^m, crs), 
derived from the joint analysis (clusters, BAO, CMB, 
Hubble parameter, and SNe) are: 

wo = -l.Utl-\t (la C.L.) , (72) 

Wa - 0.85l°;^g (1(7 C.L.) , (73) 

r!„, ^ 0.28^°;°^ {la C.L.) , (74) 

as = 0.73l°;°i^ (1^ C.L.) . (75) 

The joint analysis is compatible (at 2cr C.L.) with a cos- 
mological constant as dark energy, namely {wo,Wa) — 
(—1,0), and in agreement with constraints on Q,m and 
(Tg found by using different data on galaxy clusters (see 
discussion above) and/or different cosmological probes. 



Moreover, we find 

Fo = 69.lil;^ km/s/Mpc (la C.L.) , (76) 



in agreement with a recent determination of the Hubble 
constant from the Hubble Space Telescope Hq — 73.8 ± 
2.4 km/s/Mpc [54], and C = -0.08tJ]:J5g (la C.L.). 

The results of the global analysis are practically inde- 
pendent of n. 

Figure 4 shows the la, 2a, and 3a confidence level 
contours in the {wo,Wa) plane superimposed on the re- 
gions representing different types of evolving dark en- 
ergy models according to the classification of Ref. [l3| 
(see discussion in Section II). At the la level, phan- 
tom models of evolving dark energy of top-down type 
are slightly favored over cooling models and considerably 
preferred over thawing, pure, and bottom-up phantom 
models. Non-phantom barotropic fluids are ruled out. 
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FIG. 5: la, 2a, and 3a confidence level contours in the (fim, w) plane (left panel) and in the {Qrn,as) plane (right panel), for 
the XCDM parametrization, determined from galaxy cluster number count (including systematic errors), BAO, CMB, Hubble 
parameter, and SNe observations. Filled contours are the result of combining BAO, CMB, Hubble parameter, and SN data, 
empty (open) thin contours give the confidence level contours for cluster (CL) data only, while empty (closed) thick contours 
come from the combination of all data. 



Va. Case w ~ const 

Figure 5 shows the results of the analysis for the case 
Wo — w — const and Wa = (the XCDM parametriza- 
tion), in the {fl„i,w) and (J7,„,cr8) planes. 

Only clusters- Cluster data alone put very weak 
bounds on w and, as in the case of general evolving dark 
energy with parameters {'Wo,'Wa), prefer large values of 
rim I and relatively small values of a^'- 

n,n > 0.43 (Ict C.L.) , (77) 
as < 0.67 (Icr C.L.) . (78) 

The right panel of Fig. 5 shows the correlation between 
flm and fjg that can be approximatively parameterized 
by 

agn^J^ = 0.49 ± 0.05 {2a C.L.) , (79) 

with a slightly smaller error on a^ with respect to 
Eq. JSSl). 

Combined data analysis.- The joint data analysis gives 

w = -0.96+°;°^ (Icr C.L.) , (80) 

n„, = 0.28^°:°^ (1(7 C.L.) , (81) 

as = 0.73+°;°i^ (Icr C.L.) . (82) 

Also, we find 

Ho = 69.0li;4 km/s/Mpc (Icr C.L.) , (83) 



and ^ = -0.08lo;go (Icr C.L.). These results are almost 
independent of n. 

Observational constraints on the XCDM parametriza- 
tion have been derived from many different data sets, 
hence it provides a useful basis for comparing the dis- 
criminative power of different data. It is well known that 
SNela apparent magnitude versus redshift, BAO peak 
length scale, and CMB anisotropy data generally provide 
the most restrictive constraints on cosniological param- 
eters. Clearly, currently available massive cluster evolu- 
tion data is nowhere near as constraining as these data. 
However, cluster data results in constraints that are com- 
parable to those that follow from angular size versus red- 
shift data [ill and lookback time data js^l, but are not 
as restrictive as those from galaxy cluster gas mass frac- 
tion measurements [58| or gamma-ray burst luminosity 
observations [59|. Over all, these constraints are approx- 
imately compatible with each other and with the ACDM 
model, lending support to the belief that we are converg- 
ing on a standard cosmological model. 

VI. CONCLUSIONS 

We have constrained dark energy and other cosmolog- 
ical parameters by using X-ray temperature data of mas- 
sive galaxy clusters in the redshift range 0.05 ^ z < 0.83 
with masses within a comoving radius of 1.5/i~^Mpc 
greater than the fiducial value 8 x 10^'^h~^MQ. 
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In these analyses we have accounted for the depen- 
dence of quantities related to cluster physics — such 
as the critical density contrast, the growth factor, the 
mass conversion factor, the virial overdensity, and the 
virial radius — on the cosmological model parameters 
WQ,Wa,^m,crs,n, and i?o- 

We have found that currently available cluster data do 
not tightly constrain the dark energy equation of state, 
and that a cosmological constant is consistent with these 
observations. The analyses show that cluster data prefer 
small values of the amplitude of mass fluctuations erg, 

only clusters: as < 0.69 (IctC.L.), (84) 

as well as larg values of nonrelativistic matter energy 
density, 

only clusters: fi„i > 0.38 (IcrC.L.) . (85) 

These bounds are in slight tension at Icr with those ob- 
tained from the 7-year WMAP observations as — 
0.80±0.03 (68%C.L.) andr2,n = 0.265±0.011 (68%C.L.) 
(these results are for a spatially flat universe with cos- 
mological constant and Hq — 71 km/s/Mpc), although 
compatible at 2a confidence level. In addition, we have 
found the following normalization of a^: 

only clusters: dg^^,!/^ = 0.49 ± 0.06 (2cr C.L.) . (86) 

Cluster data alone are not sensitive to the value of the 
index n of the power-law power spectrum of the density 
perturbations, and are only very weakly dependent on 
the Hubble constant Hq. 

In order to break the ilrn.-o's degeneracy and put 
more stringent constraints on cosmological parameters, 
we have combined cluster data with BAO peak length 
scale observations, CMB anisotropy data, Hubble con- 
stant and Hubble parameter measurements, and type la 
supernova magnitude-redshift observations. In this case 
we find 

all data: as ^ 0.73tom (IcrCL.) (87) 

and 

all data: fl^ = O-^Sto^ol (l^CL.), (88) 

which are in good agreement with previous constraints 
in the literature (such as those coming from WMAP). 

Regarding the equation-of-state parameters of dark en- 
ergy, we find 

all data: wq = -l.UtoH , Wa = O.SSlJ^'^g {la C.L.) , 

(89) 

which indicates that the joint analysis is consistent with 
a cosmological constant. Moreover, the combination of 
all data is almost insensitive to n, and constrains the 
Hubble parameter to the range, 

all data: i/o = 69.ll} j km/s/Mpc (Icr C.L.) , (90) 



consistent with recent bounds from Hubble Space Tele- 
scope observations. 

Similar results are found in the case of constant 
equation-of-state parameter time-varying dark energy 
(the XCDM parametrization) . 

Our results suggest that, among models of dark en- 
ergy with varying equation of state (i.e., with Wa ^ 0), 
the top-down phantom models, for which the equation of 
state crossed the phantom divide line from a higher value 
to a lower value, are preferred over non-phantom thaw- 
ing and cooling models. Finally, non-phantom barotropic 
fiuids are excluded as models of dark energy. 

While currently available massive cluster data do not 
constrain cosmological parameters as tightly as do SNela 
apparent magnitude versus redshift measurements, or 
CMB anisotropy data, or BAO peak length scale obser- 
vations, the cluster measurements do provide constraints 
comparable to those from some of the other available 
data sets. They also play a useful role in constraining 
cosmological parameters when used in conjunction with 
other data. More importantly, we look forward to su- 
perior quality near-future massive cluster data, and an- 
ticipate the significantly more restrictive parameter con- 
straints that will result from using the techniques we have 
developed here. 
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APPENDIX A. CRITICAL DENSITY 
CONTRAST AND GROWTH FACTOR 

Critical Density Contrast.- The critical density con- 
trast Sc depends on the redshift and on the cosmology 
and can be evaluated, using the approach of Ref. [60| , as 
follows. Consider the full nonlinear equation describing 
the evolution of the density contrast: 

where a prime indicates differentiation with respect to 
the scale factor a and 

^ ^Pm Pclustcr Pm (92) 

Pm Pra 

is the density contrast with pciustor being the cluster mat- 
ter density. Since the above equation describes the non- 
linear growth of the density contrast, its value at some 
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chosen collapse time tcoiiapsc diverges. The critical den- 
sity contrast Sc at the time ^collapse is, by definition, the 
value of the density contrast at the time icoiiapsc obtained 
by solving the linearized version of Eq. (j9ip , namely 



3 3 flm 



(5 = 



(93) 



with boundary conditions for S such that, the same 
boundary conditions applied to the nonlinear equa- 
tion (I9ip makes S divergent at icoiiapsc- Following 
Rcf. 60], wc take for the initial derivative of 5, S'{ai), 
the value S'{ai) — 5 x 10~^ where = 5 x 10^^, while 
the initial value of the density contrast, S{ai), is found 
by searching for the value of S{ai) such that S diverges 
at the time tcoiiapso- We assume (as in Ref. [g^) that the 
divergency is achieved, numerically, when S exceeds the 
value 10''. 

In Fig. 1, we plot the critical density contrast, dc, as 
a function of the redshift for different values of {wq, Wa)- 



For large redshifts — where the effects of dark energy be- 
come negligible compared to those of nonrelativistic mat- 
ter — the Universe effectively approaches the Einstein-de 
Sitter model where the critical density contrast is inde- 
pendent of the redshift and is 5c = (3/20)(127r)2/3 ~ 
1.686 [sj. In Fig. 6, we show the ^c-isocontours in the 
(wo, Wa)-plaiie for different values of the redshift and for 
= 0.27. 

Growth Factor - The z-dependent part of the matter 
power spectrum — the growth factor D[z) — is 



D{z) 



5{Q) 



(94) 



and satisfies the linearized equation (j93p . The boundary 
conditions we impose are D{a = 1) = 1 and D{a = a^) = 
ai, where as before = 5 x 10~^. In Fig. 1 we plot 
the growth factor, D{z), as a function of the redshift for 
different values of {wo,Wa)- 
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APPENDIX B. MASS CONVERSION 

The virial mass in the PS or ST parametrization needs 
to be expressed as a function of the observed mass M' 
within a reference comoving radius of i?Q = 1.5/i~"'^Mpc. 
As described in Section Illb, what we really need is 
the fiducial virial mass Mq as a function of the fidu- 
cial mass within i?g adopted in the observation, namely 
M() = 8x 10"/i-^Mq. We use the following procedure 
to accomplish this. We first determine the physical virial 
radius r^^o, within which the virial mass Mq is contained, 
through the relation 

Mo{z) = —rl,iz)p„,{z)A,{z) ^ — Ri.pl^^ A,{z) , 

(95) 

where R^fi = {1 + z) Vy^, and A„ is the virial overdensity 
and is discussed in Appendix C. We then scale the virial 
mass Mq to the 1.5ft-^^Mpc comoving radius assuming a 
Navarro-Frenk- White profile for the virialized halo mass 
density [6l|: 



Pclustcr(f) 



'cluster 



(r/r,)(l + r/r,)2 



(96) 



where r is the physical radial distance and rg is a physical 
scale radius. 

Technically the procedure is as follows. From Eq. ([M)) 
we can obtain the mass M(< r) contained in the physical 
radius r: 



M{< r)^4TT / dr'r'2pciustcr(r') 
= l6TTpciustcr{rs)r^f{rs/r) 



where 



In 



1 + a; 



1 + x 



(97) 



(98) 



Both Ry,o and c„_o in Eq. (1991) are functions of Mq. 
From Eq. (|95|) we can express Ry^ as a function of Mq. 
For Cy^o, we proceed as follows. First we consider the ex- 
pression for the concentration parameter found by Bul- 
lock et al. [ill in their A-body simulation: 



CyiMy) = 



B 



l + z 



My{0) 



1 -0 



M,(0) 



(102) 



Here B ~ 9, /3 ~ 0.13, ^'^ My is the mass within a physical 
virial radius r^, 



My{z) 



An 



rl{z)prniz)Ay{z) 



(103) 



and M, (z) is a fiducial mass defined by 

a{R4z),z)^6ciz) , (104) 

where the comoving radius R*{z) is defined by 

47r 



M^iz) = ^ Rl{z) p^^ 



(105) 



Then, evaluating Eq. (jl02p for My = Mq [which corre- 
sponds to evaluating Eq. (jlOip for r„ = ry ^], we have 



where 



Cyfi = Ci,(A/o) 



B 



l + z 



Mo 
K 



M'q 



1 -0 



M,(0) 



(106) 



(107) 



is, formally, the concentration parameter (|102l) evaluated 
at [see Eq. (fT02)) ]. and 



7(z) 



Ay{z) 



1-0 



Ay{0) 



(108) 



Then applying Eq. (|97l) to the mass Mq within the phys- 
ical virial radius Vy^Q, and to the mass Mq within the 
physical radius = Rq/{1 + z) corresponding to the 
comoving radius i?Q, and taking the ratio of these two 
equations, we get 



Mo 



^11,0 



/(l/c..o) 

fiRvfi/CyfiRo) 



Here 



Cv,0 



is the "concentration parameter" 



Cy 



(99) 



(100) 



(101) 



Taking into account Eqs. ^U^, and pTOl) . we find 

that the quantity A/*(0) does not depend on Mq and is 
defined by 



I](i?40)/i?o) 



0-8 



(109) 



Finally, inserting Eqs. ^ and p06)) in Eq. ((991), we 
obtain the equation that gives Mq as a function of Mq 
[namely the equation defining the function g in Eq. ([39]) ]: 



/ Ri. 



7 fMQ 



v,0 \K 



l3+l/3\ 



Rhfi 



Mq 

VI 



(110) 



— the ratio between the physical virial radius ry and the 
physical scale radius Ts — evaluated at the physical virial 
radius r„ o ■ 



To be precise, B and j3 can (weakly) depend on the cosmology 
and the values used here are those found in Ref. [62l | where a 
cosmology with cosmological constant was assumed. 
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where 



i?l.5(^) 



1/3 



1.5/i-iMpc 



(111) 



is the comoving virial radius corresponding to the mass 
Mq, normahzed to 1.5/i~^Mpc. 

In a fiat universe dominated by nonrelativistic matter 
(the Einstein-de Sitter model), the comoving virial ra- 
dius corresponding to the mass Mq — 8 x 10^^/i~^Mq is 
about 1.57ft.^^Mpc, which gives i?i.5 ~ 1. This in turns 
means, using Eq. (jllOl) . that Mq/Mq ~ 1, namely the 
virial mass contained in a sphere of comoving virial radius 
of 1.5/i~^Mpc is approximatively Mq. Due to the pres- 
ence of the third root in Eq. (Illip . the function i?i.5(z) 
is fairly insensitive to changes in the adopted cosmology 
[through and A„(z)]. Therefore, also in the case of a 
generic cosmology with evolving dark energy, we expect 
values of Mq near to Mq. 

In Fig. 1, we plot the ratio Mq/Mq as a function of the 
redshift for different values of (wq^Wo). Mq, as argued, 
turns to be of order (of a few times) Mq. 



APPENDIX C. VIRIAL OVERDENSITY AND 
VIRIAL RADIUS 

Virial Overdensity.- The virial overdensity is the ratio 
of the cluster mass density and the background matter 
density at the time of virialization: 



Ay{z = Zy) = 



Pclustcr('2i;) 
Prn (-2^1! ) 



(112) 



where Zy is the redshift at the time of virialization. 
Using the fact that the cluster mass density is 



Pclustcr(z) 



3Af, 



cluster 



A7rr^{z) 



(113) 



where Mdustcr is the halo cluster mass and r{z) the phys- 
ical halo radius, we can rewrite Eq. (|112l) as 



(114) 



where rta and are, respectively, the physical radii of 
the halo cluster at turn-around and at virialization (r^ 
is, in other words, the physical virial radius of the halo 
cluster), while Zta and 



Pclustcr (^ta) 
Pmizta) 



(115) 



are, respectively, the redshift and virial overdensity at 
the time of turn around. 

The redshift at the virialization time is 



— -^-collapse; 



(116) 



and follows from the standard assumption that clusters 
virialize at the time of collapse. 

To find the redshift at the time of turn-around, zta, 
and to obtain the virial overdensity at the turn-around 
time, we follow the procedure of Ref. [g^. First, they 
observed that the quantity {S + l)/a'^, where 6 is the 
density contrast that satisfies the nonlinear Eq. is 
proportional to l/r^{z): 



1 



3M 



cluster 



1 



47rp. 



(0) 



\z) 



(117) 



where r{z) is the collapsing sphere's radius. [It is 
straightforward to get the above equation by using 
Eqs. and (|113|) .] Then since r{z) assumes the max- 
imum value at the turn-around time, the time of turn- 
around is found by minimizing the quantity {S + l)/a^, 
where 5 is the solution of Eq. (|9ip obtained by imposing 
the boundary conditions used in Appendix A. 

Once the turn-around time is found, the virial over- 
density at the turn-around time is simply given, from 
Eqs. (1121) and ([TT5|) . by 



C - <5(zta) + 1 , 



(118) 



where S is the solution of Eq. (|9T|) . 

The ratio of the virial radius to the turn-around ra- 
dius, Tu/rta, as a function of cosmological parameters is 
analyzed below. 

In Fig. 1 we plot the virial overdensity at the time 
of virialization. Ay, as a function of the redshift for 
different values of {wQ,Wa). In an Einstcin-de Sitter 
model (i.e., for D,m = 1), the standard assumption that 
t{zy) ~ 2t{ztn) together with the fact that — rta/2 and 
( = (37r/4)2, gives A„ ~ IStt^ ~ 177.653 (independent of 
the redshift). Indeed, each curve in Fig. 1 approaches this 
value for a sufficiently large value of the redshift since, 
as already noted, the Universe then enters the Einstein- 
de Sitter regime where the effects of dark energy become 
subdominant with respect to those of nonrelativistic mat- 
ter. 

In Fig. 7, we show the A„-isocontours in the {wo,Wa)- 
plane for different values of the redshift and for ftm — 
0.27. 

Virial Radius.- In order to "estimate" the quantity 
Tv/fta., we apply energy conservation and the virial the- 
orem to the spherical collapse of the cluster halo. 

We start by considering the total gravitational poten- 
tial energy U{r) of a sphere of radius r containing the 
cluster mass Mdustcr and dark energy: 



where 



U{r,z) = Ummir) + Umi)Eir,z) 



5 r 



(119) 



(120) 
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is the familiar gravitational potential self-energy of a 
sphere of nonrelativistic matter, and 



Ur, 



DE 



(r, z) 



47r 



[l+3w{z)]p^E{z)r^ (121) 



is the gravitational potential energy of interaction be- 
tween nonrelativistic matter and dark energy. 

Since the potential energy UmUE{r, z) depends explic- 
itly on the time, the system under consideration is not 
conservative. Therefore, neither energy conservation nor 
the virial theorem can be applied. 



^'^ Here, we are considering the case where the dark energy does 
not cluster and does not virialize, so that the only terms in the 
potential energy which are relevant for energy conservation and 
virialization are those in Eq. I)119| l [63l | . 

It is straightforward to show that the only case where UmUE (^i ^) 



In order to get a conservative system, Wang [63| 
has suggested replacing the z-dependent quantity [1 -I- 
3w(z)] pbe{z) with the same quantity evaluated at the 
turn-around time. 

Here we propose defining an effective potential energy, 
which does not depend explicitly on the time, as 



C/(<=ff)(r)={/„™(r) + t/i^«j,(r) , 



(122) 



where 

.(cff 



An 



GMcluster ( [1 + M^)] PBEiz) ) 



(123) 

is the effective potential energy of interaction and (...) is 
an operator that when applied to a z-dependent function 



is 2-independent, and the system is conservative, is that of the 
cosmological constant, namely w(z) = —1 for all times. 
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Tpiz) gives a z-independent quantity, i.e. 



dz 



= 



(124) 



The action of the (...) -operator is specified below. 

The introduction of the effective energy potentials 
()122p and (|123p . allow us to use the energy conserva- 
tion theorem that, applied at the times of virialization 
and turn-around, gives 



/C(r„) + C/(°ff)(r,) = C/(°ff)(rta) 



(125) 



where K,{r^) is the kinetic energy at the virialization time. 
Using the virial theorem 



r rfC/(<=ff)(r) 



(126) 



2 dr 

the energy conservation equation (jl25p takes the form 



(127) 

Taking into account Eqs. (|113p and (|115p . the mass of 
the cluster is 



-^^clustor - 

so that Eq. (fm]) reads 



(128) 



4??g(l+3u;ta) — -2 [1 + 775(1 + 3«;ta)] — +1 = 0. 

V na / V rta. J 

(129) 

Here Wta = w{zta) and we have defined, following 
Ref. IgI], the quantity 



CPrn{zta) 



(130) 



We have also introduced the "deviation parameter", rj, 

as 



{[l + 3wiz)]puE{z)) 

(1 + 3wta) /5DE(zta) 



(131) 



It is worth noting that in the case of a cosmological 
constant the system is conservative (see footnote 15), 
and the equation determining r^/rta is formally given 
by Eq. (|129p with 77 = 1. Hence, the only restriction to 
the action of the (...)-operator is that, when applied to 
the function [1 + 3w{z)] py)e{z), it must give 77 = 1 for 
w{z) = —1. 

Taking {i]j{z)) = tp{zta) corresponds to the choice of 
Wang, which also implies 77 = 1 for equation-of-state pa- 
rameter w{z). However, taking {ip(z)) = 'ip{z), where z 
can be anywhere in the interval [zy, zta], is also a plausi- 
ble choice. 

In the upper panel of Fig. 8, we plot the ratio of the 
virial radius to the turn-around radius in the case 77 = 1 




FIG. 8: Upper panel. The ratio of the virial radius to the 
turn-around radius, x = r„/rta with rj = 1 (corresponding to 
the Wang's choice), as a function of the redshift for different 
values of {wo,Wa) (the same as in Fig. 1). Lower panel. The 
ratio of X to y, where y = r„/rta with r] defined in Eq. p32[) . 
In both panels, we fixed Sim = 0.27. 



(Wang's choice) as a function of the redshift for different 
values of (twq, Wa). In the lower panel, we show the same 
ratio for the case ( ip{z) ) = tp(zy). For notational clarity, 
we indicate those ratios by 



rv_ 

rta 



if 77 = 1 



y if »y 



(1 + 3mi„) p-DEjZv) 
(l+3-!i)ta) PDE(zta) ' 



(132) 



where Wy = w{zy). For z 3> 1, where dark energy effects 
can be neglected, the ratio of the virial to the turn-around 
radii approaches, as expected, the asymptotic Einstein-de 
Sitter value x — 1/2. The variations of x with the choice 
of the functional form of the deviation parameter 77 are 
of order of a few percent. The resulting analysis on the 
growth of massive galaxy clusters does not appreciably 
depend on 77, and the results presented in the previous 
sections are for the case rj = 1. 
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